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Abstract. This paper presents a method to determine effective temperatures, angular semi-diameters and bolometric corrections 
for population I and II FGK type stars based on V and 2MASS IR photometry. Accurate calibration is accomplished by using 
a sample of solar analogues, whose average temperature is assumed to be equal to the solar effective temperature of 5777 K. 
By taking into account all possible sources of error we estimate associated uncertainties better than 1 % in effective temperature 
and in the range 1.0-2.5% in angular semi-diameter for unreddened stars. Comparison of our new temperatures with other 
determinations extracted from the literature indicates, in general, remarkably good agreement. These results suggest that the 
effective temperaure scale of FGK stars is currently established with an accuracy better than 0.5%— 1%. The application of the 
method to a sample of 10999 dwarfs in the Hipparcos catalogue allows us to define temperature and bolometric correction (K 
band) calibrations as a function of (V - K), [m/H] and logg . Bolometric corrections in the V and K bands as a function of T c g, 
[m/H] and logg are also given. We provide effective temperatures, angular semi-diameters, radii and bolometric corrections in 
the V and K bands for the 10999 FGK stars in our sample with the corresponding uncertainties. 
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1. Introduction 

Effective temperature and luminosity are two fundamental stel- 
lar parameters that are crucial to carry out tests of theoretical 
models of stellar structure and evolution by comparing them 
with observations. The accuracy in the determination of other 
stellar properties, such as metallicity, age or radius, hinges on 
our ability to estimate the effective temperatures and luminosi- 
ties. 

There are several approaches in the literature to compute 
effective temperature and/or luminosity. Except when applied 
to the Sun, very few of them are direct methods that permit 
an empirical measurement of these parameters. Usually, semi- 
empirical or indirect methods are based to a certain extent 
on stellar atmosphere models. Among the direct approaches 
we find the remarkable work by Code et al. (1976), which is 
based on interferometric measurements of stellar angular semi- 
diameters (9) and total fluxes (F\, {) at Earth, and the more re- 
cent works of Mozurkewich et al. (2003) and Kervella et al. 
(2004), also based on interferometry. On the other hand, indi- 



rect methods are mainly based on the use of photometry, spec- 
troscopy, or a combination of both. In the case of the temper- 
atures, although many of the published calibrations claim to 
have uncertainties of the order of several tens of degrees, val- 
ues obtained by different authors can easily have discrepancies 
of 100 K or even larger. The reason for such nagging differ- 
ences must be found somewhere in the ingredients of the meth- 
ods: atmosphere models, absolute flux calibrations, oscillator 
strengths, calibration stars, etc. 

In this paper we present a semi-empirical method to deter- 
mine effective temperatures (T e ff) and bolometric corrections 
(BC) from 2MASS 1 JHK 2 photometry (Cutri et al. 2003) that 
is applicable to FGK type stars. As all others, our method is 
susceptible to problems derived from the uncertainties in the 
ingredients mentioned above. However, our approach benefits 
from two major features: First, it provides a way to evaluate 
realistic individual uncertainties in T e ff, 6 and luminosity by 
considering all the involved errors; and second, as it is cali- 
brated to use the 2MASS photometry, it allows the calculation 
of consistent and homogeneous T e g and BC for several million 
stars in the 2MASS catalogue. This paper also provides r e if, 



* Table 3 is only available in electronic form at the CDS via 
anonymousftptocdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat?] /A+A/ 



http : //www . ipac . caltech.edu/2mass 
Throughout the paper, K refers to K s band. 
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angular semi-diameters, radii and BCs for 10999 dwarfs and 
subdwarfs in the Hipparcos catalogue ESA (1997). Such large 
sample has allowed us to construct simple parametric calibra- 
tions as a function of (V - K) , [m/H] and logg. Note that a 
preliminary version of the method presented here was already 
successfully applied to the characterization of the properties of 
planet-hosting stars (Ribas et al. 2003). 

The present paper is organized as follows. Section 2 
presents the method and explains in detail the procedure to ob- 
tain r e ff and angular semi-diameters, including the fitting al- 
gorithm, zero point corrections and error estimates. The com- 
parison of our temperatures with several previous works, both 
based on photometric and spectroscopic techniques, is de- 
scribed in Sect. 3. In Sect. 4 we present simple parametric 
calibrations of r e tf and BC as a function of (V - K)o, [m/H] 
and logg valid for dwarf and subdwarf stars. The sample of 
10999 stars used to build the calibrations is also described in 
this section together with a detailed explanation of the different 
contributors to the final uncertainties. Finally, the results are 
discussed in Sect. 5 and the conclusions of the present work 
are presented in Sect. 6. 

2. The Spectral Energy Distribution Fit (SEDF) 
method 

The use of infrared (IR) photometry to determine effective tem- 
peratures was initially proposed by Blackwell & Shallis (1977). 
Their so-called Infrared Flux Method (IRFM) uses the ratio be- 
tween the bolometric flux of the star and the monochromatic 
flux at a given infrared wavelength, both measured at Earth, 
as the observable quantity. This ratio is then compared with 
a theoretical estimate derived from stellar atmosphere models 
to carry out the determination of the effective temperature. The 
IRFM has been widely used by a number of authors, being most 
noteworthy the work by Alonso et al. (1995, 1996a,b). 

The Spectral Energy Distribution Fit (SEDF) method that 
we propose here follows a somewhat different approach, 
namely the fit of the stellar spectral energy distribution from the 
optical (V) to the IR (JHK) using synthetic photometry com- 
puted from stellar atmosphere models. Unlike the Alonso et al. 
(1996a) implementation of the IRFM, which averages temper- 
atures derived individually for each IR band, our method takes 
into account the four bands simultaneously (and naturally). 
In addition, and also unlike the IRFM, the bolometric flux is 
not required a priori by the SEDF method but results self- 
consistently with the temperature. The fitting algorithm (see 
Sect. 2.2) minimizes the difference between observed and syn- 
thetic photometry by tuning the values of the effective temper- 
ature and the angular semi-diameter. The BC can be obtained 
from these two parameters, and then, when the distance to the 
star is known, the luminosity is computed from the BC and the 
absolute magnitude in a given photometric band. The uncer- 
tainties of the derived parameters (T e ff , angular semi-diameter 
and BC) are estimated from the errors in the observed and syn- 
thetic photometry as well as in the assumed [m/H], logg and 

Ay. 

From a theoretical point of view, the SEDF method could 
be applied to stars of any spectral type and luminosity class. 



However, the IR flux becomes very sensitive to metallicity and 
surface gravity for stars hotter than 8000 K so that small uncer- 
tainties in these parameters translate into large uncertainties in 
the effective temperature. In such situation the SEDF approach 
becomes inadequate. At the cold end, the accuracy of stellar 
atmosphere models limits the use of the method to stars hotter 
than 4000 K (molecular opacity plays an important role below 
this temperature). These limitations restrict the applicability of 
the SEDF to FGK type stars. Fortunately, these stars are very 
common in the Galaxy and dominate the content of most of 
survey catalogues. They are crucial for several key astrophys- 
ical topics, such as the study of the structure and evolution of 
the Galaxy, both the disk and the halo, and the characterization 
of planet-hosting stars, among others. 

2.1. Calculation of synthetic photometry 

The calculation of the synthetic photometry requires a well- 
characterized photometric system, an accurate flux calibra- 
tion and suitable synthetic spectra. The work by Cohen et al. 
(2003a,b) provides consistent absolute flux calibrations in both 
the visible (V) (Landolt system) and IR (2MASS JHK) bands. 
The calibration given by Cohen et al. is computed from a set 
of calibrated templates, using the synthetic Kurucz spectrum 
of Vega of Cohen et al. (1992). In the case of the IR photom- 
etry, they consider the transmission of the camera and filters, 
the detector properties and the Earth's atmosphere character- 
istics. From the comparison between observed and synthetic 
photometry for a set of 9 A-type stars and 24 cool giants, the 
authors infer the need to introduce a zero point offset in the syn- 
thetic photometry to match the observed 2MASS photometry: 
0.001+0.005 mag (/); -0.019+0.007 mag (H); 0.017+0.005 
mag (K). The calculation of such values is not exempt of some 
difficulty since the dispersions of the differences between both 
photometries (synthetic and observed) are of the same magni- 
tude as the zero point offsets themselves. 

To compute the syntheric magnitudes we made use of the 
no-overshoot Kurucz atmosphere models grid (Kurucz 1979) 
taken from http : //kurucz . harvard . edu/grids . html: 



KynCTrff, logg, [m/H]) = 2.5 log 



F' 



{Fi(T eS ,logg, [m/H])) 



(1) 



where F' caJ is the absolute flux calibration given by Cohen et al. 
(2003b) (for = 0) and Fi(T eB ,g, [m/H]) is the flux in the 
;' band computed from the integration of the model atmosphere 
convolved with the transmission function (filter, detector and 
Earth's atmosphere) from Cohen et al. (2003b): 



FiiT^Xogg, [m/H]) 



Jo 



0(r eff , logg, [m/H], A)Ti(A)dA (2) 



where (f>(T e g, logg, [M/H], A) is the flux given by the stellar 
atmosphere model and 71(A) the effective transmission function 
in the ;' band normalized to a peak value of unity. 

2.2. Fitting algorithm 

The fitting algorithm is based on the minimization of the x 1 
function defined from the differences between observed (cor- 
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rected for interstellar extinction) and synthetic VJHK magni- 
tudes, weighted with the corresponding error: 

^ = | V-Ay- V syn j 2 + ^-A./-y syn j 2 + 

+ (H -Ah - H syn \ 2 ^ / K - A K - K sya \ 2 

I o-h ) \ o-k ) 

This function depends (via the synthetic photometry) on T e ff, 
\ogg, [m/H] and a magnitude difference 3K, which is the ratio 
between the synthetic (star's surface) and the observed flux (at 
Earth) ( J[ = -2.5 \og F star /F Ea rth)- & is directly related to the 
angular semi-diameter by the following expression: 

= iq-ow (4) 

Although the synthetic photometry depends implicitly on grav- 
ity and metallicity, in practice, the spectral energy distribution 
in the optical/IR for our range of temperatures is only weakly 
dependent on these quantities. This fact makes it possible to 
obtain accurate temperatures even for stars with poor determi- 
nations of log g and [m/H]. 

As can be seen, the x 2 function depends also on the inter- 
stellar absorption Ay (the absorption in the other bands can be 
computed using the extinction law of Schaifers & Voigt (1982): 
Aj = 0.30Ay, A H = Q.2AA V and A K = 0.15A V ). In principle, 
it is possible to consider Ay as a free parameter. However, the 
strong correlation between r e tf and Ay, especially for the hotter 
stars, decreases the precision in the determination of both pa- 
rameters, with resulting typical uncertainties of 4% in T e g and 
0.25 mag in Ay. Thus, for best performance, A v should only be 
considered as a free parameter when its value is suspected large 
and no other method for its estimation is available. In general, 
the best approach is to fix the value of A v in Eq. (3) from the 
estimation by photometric calibrations, for instance. 

Therefore, the only two adjustable parameters by the SEDF 
method in the present work are T e s and J{, whereas logg, 
[m/H] and A v are fixed parameters. To minimize Eq. (3) 
with respect to these two parameters we use the Levenberg- 
Marquardt algorithm (Press et al. 1992), which is designed to 
fit a set of data to a non-linear model. In all our tests, conver- 
gence towards the minimum value of x 2 was reached rapidly 
and unequivocally. 

2.3. Calibration of the SEDF method using solar 
analogues 

The standard procedure for the calibration of an indirect 
method to determine effective temperatures is based on the 
comparison of the results with accurate temperatures from di- 
rect methods for a set of stars. In this way, the list of stars with 
empirical effective temperatures and angular semi-diameters 
from Code et al. (1976) has been widely used for calibration 
purposes. This list has been recently increased with the the 
works of Mozurkewich et al. (2003) and Kervella et al. (2004). 
Other authors use well-studied stars, such as the Sun, Vega or 
Arcturus, to calibrate their methods. 

Unfortunately, the few stars with empirical values of 
r e ff are too bright to have accurate 2MASS photometry and 



they are of no use to calibrate the SEDF method. As an alter- 
native, we have used the list of photometric solar analogues 
compiled by Cayrel de Strobel (1996). We assume that, as an 
ensemble, the average of the effective temperatures of these 
photometric solar analogues should be equal to the effective 
temperature of the Sun (i.e., 5777 K). 

After selecting a subsample of 50 unreddened stars 
with non-saturated 2MASS photometry from table 1 of 
Cayrel de Strobel (1996), we computed their temperatures us- 
ing the SEDF method. We obtained an average temperature of 
5832 ± 14 K, i.e., 55 K (or -1%) higher than the solar effec- 
tive temperature. Exactly the same value is obtained if we use 
the subset of solar "effective temperature analogues" from ta- 
ble 5 of Cayrel de Strobel. Without a profound analysis of all 
the ingredients involved - from the stellar atmosphere model to 
the absolute flux calibration, - it is very difficult to assess the 
reasons for such difference. However, it seems clear that the 
temperature scale as obtained from the synthetic photometry 
alone needs a correction to agree with the average of the solar 
analogues. From a formal point of view, this correction can be 
computed from the synthetic photometry that results from forc- 
ing a value of T e g = 5777 K to the entire sample. After doing so, 
we replaced the zero points given by Cohen et al. (2003b) (see 
Sect. 2. 1) by the average difference (for each band) between the 
observed and synthetic photometry computed for the solar ana- 
logues. Assuming that there is no offset on the V band, the off- 
sets for the other bands are: 0.027+0.003 mag (J); 0.075±0.005 
mag (H); 0.022+0.005 mag (K). It is interesting to note that 
both in the case of Cohen et al. (2003a) and in our case, the 
value of the offset in the H band differs significantly from the 
offsets in J and K. It should be stressed that the effective tem- 
peratures given by Cayrel de Strobel ( 1 996) have not been used 
here. We have only used the property of the stars in being clas- 
sified as solar analogues, and, consequently, we assumed their 
average temperature to be equal to the solar effective tempera- 
ture. 

In our procedure, we are implicitly assuming that the cor- 
rection in our temperature scale is just a zero point offset and 
that no dependence on temperature or metallicity is present. 
These assumptions are justified a posteriori in Sect. 3, where 
several comparisons of SEDF temperatures with other photo- 
metric and spectroscopic determinations are shown. 

The angular semi-diameters computed from Eq. (4) were 
used to check the consistency of the new zero points in our 
temperature scale. These angular semi-diameters were com- 
pared with the direct values compiled in the CHARM2 cata- 
logue (Richichi & Percheron 2005). We restricted the compar- 
ison to stars with accurate VLBI or indirect (spectrophotome- 
try) measurements of the semi-diameter. Only 10 of these stars 
fulfill the conditions for applicability of the SEDF method. 
Figure 1 shows the comparison of the semi-diameters for these 
10 stars. The agreement is excellent, with an average differ- 
ence (%,- - #sedfX weighted with the inverse of the square 
of the error, equal to -0.3% with a s.d. of 4.6% (see Table 
1). All the direct values used in the comparison correspond 
to an uniform stellar disk. A crude comparison of both uni- 
form disk and limb darkened values for about 1600 F, G and 
K stars in the CHARM2 catalogue indicates a ~4% positive 
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2.4. Error estimation 




0.2 0.3 0.4 

Semidiameter (SEDF) (mas) 



0.5 



Fig. 1. Comparison of angular semi-diameters computed from the 
SEDF method (with the new zero point in the temperature scale) and 
from the CHARM2 catalogue. In the case of HD 209458, the compar- 
ison of the semi-diameter is between the SEDF method and an empir- 
ical determination from a high-precision transit light curve. 



Table 1. Comparison of direct and SEDF angular semi-diameters. 
6* dir are VLBI and spectrophotometry values (Richichi & Percheron 
2005). For HD 209458, dir is derived from a planetary transit light 
curve. 



Star 


6W (mas) 


fsEDF (mas) 


A0(%) 


HIP 6702 


0.190±0.050 


0.189±0.004 


0.5 


HIP 8433 


0.225±0.050 


0.286±0.005 


-27.1 


HIP 48113 


0.405±0.005 


0.402±0.006 


0.7 


HIP 50786 


0.260±0.010 


0.265±0.004 


-1.9 


HIP 51056 


0.290±0.030 


0.306±0.005 


-5.5 


HIP 85365 


0.380±0.015 


0.410±0.036 


-7.9 


HIP 91237 


0.210±0.015 


0.251 ±0.004 


-19.5 


HIP 96895 


0.270±0.015 


0.278±0.004 


-3.0 


HIP 96901 


0.260±0.015 


0.255±0.003 


1.9 


HIP 113357 


0.365±0.010 


0.352±0.006 


3.6 


HD 209458 


0.113±0.010 


0.115±0.002 


-1.8 



correction for limb darkening, of the same order of the disper- 
sion of the relative differences shown in Table 1. In addition, 
we compared the radius of HD 209458 - obtained with the 
Hubble Space Telescope from a high precision planetary tran- 
sit light curve (Brown et al. 2001) - with our estimation from 
SEDF, obtaining very good agreement: 1.146 ± 0.050 R Q and 
1.160 ± 0.058 R & , respectively. 



One of the features of the SEDF method is that it yields indi- 
vidual uncertainties of both r e jf and 6. The total uncertainty can 
be calculated by combining the contributions from the spec- 
tral energy distribution fit and also from the uncertainties in 
the fixed parameters, i.e., \ogg, [m/H], and Ay. Assuming null 
correlation between these different (and, in principle, indepen- 
dent) sources of error, the total uncertainty of the Y k parameter 
(r e ff or J{) is the quadratic sum of the different error contribu- 
tions: 



(AY k ) 2 = (AF, SEDF ) 2 + 



dY k 



d[m/H] 



[m/H]) 2 



dY k 
dlogg 



dY k 



(crioss) + 'a47' {(TAv) 



(5) 

As mentioned above, the minimization of^- 2 is carried out with 
respect to the four magnitudes VJHK, with the parameters A v , 
[m/H] and log g being held fixed. The first term of the equation 
(AY k SEDF ) is the error in the parameter Y k coming from the fit 
to the spectral energy distribution, which is computed from the 
covariance matrix (AY k SEDF = VC«). The derivatives in the 
other three terms of the equation are determined numerically: 

dY k ^ Y k (+A[m/H]) - Y k (-A[m/H]) 
d[m/H] 2 ( ' 

in the case of the metallicity, and in an analogous way for the 
surface gravity and the interstellar absorption. 

The error in T e g is obtained directly from Eq. (5), whereas 
the error in 8 must be calculated from the error in the J[ param- 
eter: 



ere =0.2 In 10 00-^ 



(7) 



The adopted values for the errors in the magnitudes (cr m ), 
metallicity (<T[m/h]) and surface gravity (o-\ ogg ) are discussed 
in Sect. 4.1.1. 

3. Comparison with other methods 

Five samples of FGK stars with accurate determinations of 
effective temperatures (both photometric and spectroscopic) 
were selected from the literature (Alonso et al. 1996a, Ramirez 
& Melendez 2005, Fuhrmann 1998, Santos et al. 2003 and 
Edvardsson et al. 1993) to carry out a comparison with our 
results. We put special attention in correcting for the effects 
of interstellar reddening, which could lead to systematic dif- 
ferences. For the Alonso et al. and Ramirez & Melendez sam- 
ples (the most reddened), interstellar reddening was corrected 
using the values of E(B - V) given by the authors so that 
the two temperature estimations would be directly compara- 
ble. For the other three samples, composed of stars at closer 
distances, we restricted our comparisons to unreddened stars. 
Anyhow, this meant rejecting very few stars from further anal- 
ysis. Among several papers in the literature, we chose these 
five samples because they have a minimum of 25 stars with 
non-saturated 2MASS photometry and the values of [m/H] and 
log g - needed for a consistent comparison - are provided by 
the authors. 
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3.1. Methods based on IR photometry: Alonso et al. 
(1996a) and Ramirez & Melendez (2005) 

As mentioned above, the IRFM is the most popular method to 
compute effective temperatures from IR photometry. The work 
by Alonso et al. (1996a) is undoubtedly the widest applica- 
tion of the IRFM to FGK stars. The authors computed effective 
temperatures for 462 stars with known interstellar absorption, 
surface gravity and metallicity. After selecting the stars in the 
Alonso et al. sample in the range 4000 < r eff < 8000 K and 
with errors in the 2MASS photometry below 0.05 mag, we ob- 
tained effective temperatures from the SEDF method for a sub- 
set of 315 stars. The comparison between both determinations 
of r e ff is shown in Fig. 2. The average difference Ar e ff (IRFM 
- SEDF) was found to be -67 K, with a standard deviation of 
8 1 K. The dependence of this difference on the temperature is 
not significant: 7^ FM = 1.030 r e s | DF - 239 K. As shown in the 
bottom panel of Fig. 2, there is no dependence of the tempera- 
ture difference with metallicity. 

In a recent work, Ramirez & Melendez (2005) have re- 
computed the IRFM temperatures of almost all the stars in 
Alonso et al. (1996a) using updated input data. According to 
the authors, the difference between the old and new temper- 
ature scales is not significant. They also compare their effec- 
tive temperatures with some direct determinations. The authors 
conclude that there is a systematic difference of about 40 K at 
solar temperature (in the sense IRFM Alonso - their values). The 
comparison between the temperatures of Ramirez & Melendez 
and our determinations is shown in Fig. 3. For 385 stars in com- 
mon we find A7/ eff (IRFM -SEDF) equal to -58 K (cr rdr =67 
K), and T™* M = 1.061 r s | DF - 403 K. Unlike in the case 
of Alonso et al. (1996a), the dependence of AT e g in [m/H] is 
relevant (Fig. 3, bottom panel). For [m/H] < -2.0 the tem- 
peratures from Ramirez & Melendez are clearly hotter than 
our temperatures. The same trend was found by Charbonnel & 
Primas (2005) when comparing their temperatures of 32 halo 
dwarfs (-3.5 < [Fe/H] < -1.0) with the values of Ramirez & 
Melendez. 



3.2. Other methods 

Besides the IRFM, which uses IR photometry, we have also 
compared the effective temperatures obtained using the SEDF 
method with other determinations. Two of these (Fuhrmann 
1998 and Santos et al. 2003) are spectroscopic works, while 
in another case (Edvardsson et al. 1993) the temperatures are 
based on uvby -{} photometry. The chief problem in the case of 
spectroscopic determinations is that, in general, they are mostly 
applied to bright stars, which have poor 2MASS photometry 
(the 2MASS detectors saturate for stars brighter than K « 4 
mag). This fact reduces the number of stars in the Fuhrmann 
(1998) and Santos et al. (2003) samples that can be compared 
with SEDF method. 



3.2.1. Fuhrmann (1998) 

This sample is composed of about 50 F and G nearby stars, 
both main sequence and subgiants, of the Galactic disk and 
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Fig. 2. Comparison of the effective temperatures from the IRFM and 
the SEDF method for 315 stars in the sample of Alonso et al. (1996a). 
The bottom panel shows the temperature difference as a function of 
the metallicity. 



halo. Effective temperatures were determined from fits to the 
wings of the Balmer lines. Of those stars, 24 have accurate 
2MASS photometry so that reliable SEDF temperatures can be 
derived. The comparison is shown in Fig. 4. The mean aver- 
age difference Ar e ff (Fuhrmann - SEDF) is 12 K, (cr7- c|r =45 
K), with a slight dependence on the temperature: j^ hlmlinn - 
0.895 r s | DF + 618 K. No dependence was found between Ar eff 
and [m/H] (Fig. 4, bottom panel). 



3.2.2. Santos et al. (2003) 

To study the correlation between the metallicity and the prob- 
ability of a star to host a planet, Santos et al. (2003) obtained 
spectroscopic temperatures for 139 stars based on the analy- 
sis of several iron lines. Effective temperatures for a total of 
101 stars in the sample of Santos et al. can be obtained us- 
ing the SEDF method. In this case, Ar eff (Santos - SEDF) is 
28 K, with (TT elr = 68 K, and practically independent of the 
temperature: T^ ntos = 1.053 7; SFDF - 270 K (Fig. 5). There is 
no dependence of Ar e tf with [m/H] (Fig. 5, bottom panel). 
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Fig. 3. Comparison of the effective temperatures from the IRFM and 
the SEDF method for 386 stars in the sample of Ramirez & Melendez 
(2005). The bottom panel shows the temperature difference as a func- 
tion of the metallicity. 
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Fig. 4. Comparison of the effective temperatures from fits to B aimer 
lines and the SEDF method for 24 stars in common with the sample of 
Fuhrmann (1998). The bottom panel shows the temperature difference 
as a function of the metallicity. 



3.2.3. Edvardsson et al. (1993) 

The sample of Edvardsson et al. is composed by 189 nearby F 
and G type stars. In contrast with the previous two, the effective 
temperature is not derived from spectroscopy but from uvby — B 
photometry. To do so, the authors built a grid of synthetic 
photometry using the atmosphere models of Gustafsson et al. 
(1975) and further improved it by adding several new atomic 
and molecular lines. Effective temperatures for 115 stars in 
their sample could be derived using the SEDF method. The 
average difference Ar e ff (Edvardsson - SEDF) is 10 K, with 
a dispersion of 70 K and no dependence on the temperature: 
^Edvardsson = { 006 ^sedf _ 2 7 K (Fig. 6). As in the two previ- 
ous cases, the bottom panel of Fig. 6 shows that the temperature 
difference is not correlated with [m/H]. 

4. Parametric calibrations 

The practical use of the SEDF method as it has been described 
in Sect. 2 is not straightforward since it requires the calcula- 
tion of synthetic photometry from stellar atmosphere models 
and then use a numerical algorithm to minimize the x 2 func- 
tion. Parametric calibrations (as a function of one or more pa- 
rameters) may offer a suitable means to estimate reliable ef- 
fective temperatures in cases where simplicity and speed are 



to be preferred over the best possible accuracy. In this section 
we present calibrations for both T e g and BC as a function of 
(V - K)o, [m/H]and logg. To calculate the calibrations, the 
SEDF method was applied to a sample of stars in the Hipparcos 
catalogue, as described below. Note that these calibrations are 
subject to two limitations with respect to the full SEDF method: 
First, they are simplifications since not all the available infor- 
mation is used, and second, individual uncertainties cannot be 
determined. 

4.1. The stellar sample 

We collected a sample of FGK dwarfs and subdwarfs in the 
Hipparcos catalogue, and therefore with measured trigono- 
metric parallaxes. Their V magnitudes come mainly from the 
Hauck & Mermilliod (1998) catalogue, except for those stars 
with less than two observations, where we used the Hipparcos 
catalogue. The entire sample has complete and non-saturated 
JHK photometry in the 2MASS catalogue. The metallicity 
was extracted from the compilation of Cayrel de Strobel et al. 
(2001) or computed from uvby - B photometry - either mea- 
sured from our own observations or obtained from the Hauck 
& Mermilliod (1998) catalogue -, using a slightly revised ver- 
sion of the Schuster & Nissen (1989) calibration. The range of 
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Fig. 5. Comparison of the effective temperatures from iron line fits and 
the SEDF method for 101 stars in common with the sample of Santos 
et al. (2003). The bottom panel shows the temperature difference as a 
function of the metallicity. 

metallicities covered by the sample is -3.0 < [m/H] < 0.5. 
Values of \ogg were computed from uvby - /3 photometry 
(Masana 1994; Jordi et al. 1996). Originally, the sample was 
built to study the structure and kinematics of the disk and halo 
of the Galaxy (Masana (2004)) and a full description including 
the photometry and a complete set of physical parameters will 
be provided in a forthcoming paper (Masana et al. 2006). 

In spite of the proximity of the stars (90% of them are closer 
than 200 pc), we computed individual interstellar absorptions 
from uvby - (3 photometry and corrected the observed magni- 
tudes. As discussed below, interstellar absorption is one of the 
most important sources of uncertainty in the r eff determination. 

4.1.1. Errors 

For our sample, the errors in the magnitudes, metallicity and 
surface gravity were estimated in the following manner: 

- Errors in the VJHK magnitudes: The total error in each 
magnitude was computed as the quadratic sum of the ob- 
servational error, the error in the absolute flux calibration 
and the error in the determination of the interstellar ex- 
tinction. The first one comes from the photometric cata- 
logues. However, to prevent the underestimation of the er- 
ror in the V band, usually computed from the average of 



7000 



6500 



, 6000 





5500 


s 


300 


=1 


150 


ta 




™ *§ 




H° 







-150 


a. 






-300 


g 


200 


- 




W 


100 













-100 


£ % 




H 


-200 




5500 



6000 6500 
T ff (SEDF) (K) 



7000 



1 1 


' 1 1 1 ' 
• . ?.>'. J..}-. 


, i_ 1 i 1 






- — — r 




■*'. * . ' 




1 , 


' 






[m/H] 

Fig. 6. Comparison of the effective temperatures from photometry 
and the SEDF method for 115 stars in common with the sample of 
Edvardsson et al. (1993). The bottom panel shows the temperature 
difference as a function of the metallicity. The standard deviation for 
a single star in Edvardsson et al. (1993) is 81 K. 



a few measurements (and with no evaluation of systemat- 
ics), we have set a minimum error in V equal to 0.015 mag. 
The uncertainties in the absolute flux calibration are given 
by Cohen et al. (2003a,b) and are in the range 1.5-1.7% 
(0.016-0.019 mag), depending on the band. For those stars 
affected by interstellar reddening, the uncertainty in Ay as 
derived from photometric calibrations based on uvby - (5 
photometry (Jordi et al. 1996) is expected to be of about 
0.05 mag, or ~1.5% in r e ff. 
- Errors in [m/H] and \ogg: As mentioned above, 
[m/H] was obtained, whenever possible, from spectro- 
scopic measurements, and otherwise we used photometric 
calibrations, with assigned uncertainties of 0.10 dex and 
0.15 dex, respectively. We assigned uncertainties of 0.18 
dex to log g values determined from photometric calibra- 
tions. The effect on the final effective temperatures due to 
the uncertainties of both [m/H] and logg is very small: 
an error of 0.5 dex in [m/H] has an effect in r ei f of less 
than 0.5%, whereas the same error in logg has an effect 
that ranges between 0% and 1%, depending on the value of 
r etf and logg. 

No error was attributed to the flux in the stellar atmosphere 
models. Comparisons carried out by using other stellar atmo- 
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sphere models such as those by Castelli et al. (1997) and the 
NextGen models by Hauschildt et al. (1999) show resulting dif- 
ferences in temperature below ~ 0.3% in all cases (Ribas et al. 
2003). 

An estimation of the final errors in T e ff as function of r e jf , 
[m/H] and logg is shown in Fig. 7. As can be seen, the fi- 
nal error is almost independent of [m/H] and logg, but not of 
r e ff. Hotter stars have greater uncertainties (slightly >1% for 
r et f = 7500 K) than cooler stars (0.6% for r eff = 5000 K). 
It is important to note that, in the case of reddened stars, an 
uncertainty of 0.05 mag in Ay can double the error in T e tf com- 
pared to the values in Fig. 7. For the angular semi-diameter 
the behaviour of the errors is very similar to those of the ef- 
fective temperature, with values for unreddened stars of about 
1.0-2.5%. This means that for Hipparcos stars with good paral- 
laxes, we are able to determine the stellar radii with remarkable 
uncertainties of about 1.5-5.0%. 

Figure 8 shows the cumulative histograms of the relative 
errors in effective temperature, angular semi-diameter and ra- 
dius for the 10999 stars of the sample. As can be seen, about 
85% of the stars have determinations of T e tf better than 1.1%. 
The relative error in the angular semi-diameter is also better 
than 1.5% for about 85% of the stars. In the case of the radii, 
the main contributor to the error is the uncertainty in the paral- 
lax. Even so, 50% of the stars have radius determinations better 
than 10%, and 85% of the stars better than 25%. It should be 



kept in mind that most of the stars in our sample are essentially 
unreddened, thus yielding the best possible accuracy. 

Table 3 lists effective temperatures, angular semi- 
diameters, radii and bolometric corrections in the V and K 
(2MASS) bands with the corresponding uncertainties for the 
entire sample. Using these values, we calculated simple para- 
metric calibrations of effective temperature and bolometric cor- 
rection as described below. 



4.2. Effective temperature calibration 

Although the effective temperature for FGK type stars is 
strongly correlated with the (V - K)o index (see for instance 
Alonso et al. (1996a)), it also depends weakly on the metallicity 
and surface gravity, as we mentioned in Sect. 2. Therefore, an 
empirical calibration of r e ff should include terms in all (V-K)q, 
[m/H] and log g. Furthermore, in our case the calibrations were 
constructed separately in two {V-K)q intervals. Stars departing 
more than 3cr from the fit were rejected. The resulting expres- 
sions are: 

- 0.35 <(V-K) < 1.15 (4954 stars): 

<9 eff = 0.5961 +0.1567(V -K) + 0.0309(V - K) 2 Q + 
+ 0.009[m/#] + 0.0022[m/#] 2 
+ 0.0021(V - K) [m/H] - 0.0067 logg 
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Fig. 8. Cumulative histograms of the relative error in effective tem- 
perature, angular semi-diameter and radius for the 10999 stars in the 
sample. 



0.0028 



(8) 



- 1.15 < (V-K) < 3.0 (5820 stars): 

f3 eff = 0.5135 + 0.2687 (V-K) Q - 0.0174(V- 
+ 0.0298[m/ff]-0.0009[m/#] 2 
- 0.0184(V - K) [m/H] - 0.0028 logg 
o-8 M = 0.0026 



(9) 



where f3 eff = The standard deviation of Eqs. (8) and (9) is 
about 20 K and 25 K, respectively. As shown in Fig. 9, there 
is no residual trend as a function of (V - K)o, [m/H] or \ogg. 
Equation (8) is aplicable in the range 3.25 < logg < 4.75 and 
eq. 9 in the range 3.75 < logg < 4.75. Furthermore the calibra- 
tions are valid in the ranges of colours and metallicities of the 
sample: 

-3.0 < [m/H] < -1.5 for 1.0 < (V - K) Q < 2.9 
-1.5 < [m/H] < -0.5 for 0.5 < (V - K) < 2.9 
-0.5 < [m/H] < 0.0 for 0.4 < (V - K)o < 3.0 
0.5 < [m/H] < 0.5 for 0.35 < (V - K) Q < 2.8 (10) 

While (V - K)o is an observational quantity and [m/H] can 
be obtained from photometric and/or spectroscopic measure- 
ments, a good determination of logg is usually unavailable 
for the most of the stars. This could severely restrict the ap- 
plicability of the above calibrations. However, some photo- 
metric indexes, as the Stromgren 5c\ (see Crawford (1975) or 
Olsen (1988)), are good surface gravity indicators and, if avail- 
able, can help to estimate log g . On the other hand, catalogues 
of spectroscopic metallicities usually provide an estimation of 
the surface gravity. Finally, a crude estimation of log g can be 
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Fig. 9. Residuals of the T eS fit as function of effective temperature, 
metallicity and surface gravity. 



get from MK classification. The error in effective temperature 
caused by an error in log g will be: 



AT, 



eft' 



5040 



T 2 eS Alogg 



(11) 



where a is the coefficient of the log g terms in Eqs. (8) and (9). 
In the worst case, that of the hotter stars, Ar e fj = 85 A log g. 
Thus, even if the uncertainty in log g is as much as 0.5 dex, the 
error induced in T e s is just 40 K. 

The fits for four different metallicities and log g = 4.5 to- 
gether with the stellar sample are shown in Fig. 10. Figure 11 
shows the empirical T e g -(V - K)q relationships as a function 
of the metallicity. 

4.3. Bolometric correction calibration 

Since the SEDF method provides both effective temperature 
and angular semi-diameter, it also naturally allows for the de- 
termination of the bolometric correction in a specific band. 
From this, if the distance is known, one can compute the lu- 
minosity of the star. The bolometric correction in a given band 
is defined as the difference between the bolometric magnitude 
and the magnitude in that band: 



Bd = M bo \ - Mi = OTboi - m t 



(12) 



where m\, Q \ and m, are assumed to be corrected of interstellar 
reddening. M^x can be easily expressed as a function of the 
radius and effective temperature: 



M bol = -5 log A - 10 log + 4.74 

K& * eff © 



(13) 
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Fig. 10. r eff -(V - K) fits for four groups of stars with different 
metallicities. The empirical relationships correspond to log g =4.5 and 
[m/H] =-2.0, -1.0, -0.25 and +0.25. 
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Fig. 11. r e ff -(V - K)o relationships for logg =4.5 and four different 
metallicities. 



where fl G =6.95508 10 8 m and T efs = 5777 K. For the Sun we 
adopt V(O) = -26.75 mag and m bo i(0) = -26.83 mag, and 
therefore BC V (Q) = -0.08 (Cox 2000). 

Using the definition of the absolute magnitude at a given 
band (M x — m x + 5 log n + 5) and expressing the radius as 
function of the parallax (n) and the angular semi-diameter 
(R = 6 In), we obtain the following formula for the bolomet- 



ric correction: 



BC X = M ho] - M x 



-5 log | <K— j - 10 log — - 0.26 - m x 



(14) 



where *7C is the factor corresponding to the transformation of 
units. Once the bolometric correction for a band i is known, the 
bolometric correction for any band j can be determined from: 



BCj — (mi — mj) + Bd 



(15) 



The error in the bolometric correction can be expressed as a 
function of the uncertainties in T eS , 6 (or J{) and m,, as in Sect. 
2.4: 



5 erg 
In 10 e 



10 0r eff 



1 In 10 r eff 



In 10 r eff 

2 



+ (o- m f = 



+ (cr m f 



(16) 



The procedure described here was used to compute the 
bolometric correction in the K (2MASS) band for the stars 
in our sample. In the same way as for the effective tempera- 
ture, we calibrated BC as a function of (V - K) , [m/H] and 
log g with the following results: 

- 0.35 < (V - K) < 1.15 (4906 stars): 

BC K = 0.1275 + 0.9907(7 - K) ~ 0.0395(V - K)l + 

+0.0693[m/#] + 0.0140[m///] 2 

+0.0120(V - K) [m/H] - 0.0253 logg 
o-bc = 0.007 mag (17) 

- 1.15 <(V-K)o< 3.0 (5783 stars): 

BC K = -0.1041 + 1. 2600(7 - K) - 0A570(V - K)l + 

+0.1460[m/#] + 0.0010[m/ff] 2 

-0.0631(V - K) [m/H] - 0.0079 logg 
(Tbc = 0.005 mag 



(18) 



The range of validity of these calibrations is the same as in the 
case of the effective temperature. The bolometric correction in 
any band can be obtained from BCk via Eq. (15). 

Figure 12 shows the fits for four different metallicities, to- 
gether with the stars in the sample used to obtain the calibra- 
tions. The BCk - (V - K)q relationships as a function of the 
metallicity are shown in Fig. 13. The calibration is tabulated 
in Table 2 and compared with the calibrations by Alonso et al. 
(1995) and Flower (1996) in Fig. 14, showing good agreement. 

5. Discussion 

The procedure described in this paper yields three basic stellar 
parameters: the best-fitting effective temperature and angular 
semi-diameter and, from them, the bolometric correction. If the 
distance is known, 6 can be transformed into the true stellar 
radius. The accuracies of the parameters for the stars in our 
sample are 0.5-1.3% in r eff , 1.0-2.5% in and 0.04-0.08 mag 
for the BC. 
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Comparisons with other determinations described in Sect. 
3 show general good agreement, with differences below 0.5<x, 
except for Alonso et al. (1996a) and Ramirez & Melendez 
(2005), where the difference is about 0.8cr. The use of differ- 
ent atmosphere models and the intrinsic nature of the meth- 
ods (photometric for Edvardsson et al., Alonso et al., Ramirez 
& Melendez and ours; spectroscopic for Santos et al. and 
Fuhrmann) can explain in part the small differences. In the 
case of the IRFM, the main difference between the implemen- 
tation of both Alonso et al. (1996a) and Ramirez & Melendez 
(2005), and the SEDF method is the absolute flux calibration: 
Alonso et al. (1994) for the IRFM and Cohen et al. (2003b) for 
the SEDF. This, together with the use of different versions of 
the ATLAS9 atmosphere models, is probably the reason for the 
~60 K differential between both implementations of the IRFM 
and our determination. For Ramirez & Melendez (2005) there 
is a dependence of Ar e ff with [m/H] in such a way that the 
temperature difference (Ramirez & Melendez- SEDF) increase 
abruptly for [m/H] < 2.0. In all the other cases, the temperature 
differences are not correlated with [m/H]. 

The most important factor to explain the systematics among 
the effective temperatures computed from different methods is 
the absolute flux calibration affecting photometric determina- 



tions and inaccuracies of model atmospheres (non-LTE effects, 
3D effects, treatment on convection ,...) affecting both photo- 
metric and spectroscopic determinations. Bohlin & Gilliland 
(2004) pointed out a probable 2% overestimation of the IR flux 
in the Vega model used by Cohen et al. (2003a). A 2% shift in 
absolute flux calibration is equivalent to a difference of about 
40 K in temperature and to a zero point offset in the synthetic 
photometry of 0.022 mag. Such value would be compatible 
with our magnitude zero points in Sect. 2.3. 

Beyond the internal errors, which in the case of the SEDF 
take into account the uncertainty in the flux calibration and all 
other error sources, the comparison with other methods shows 
that, at present, the systematic errors involved in the determi- 
nation of effective temperature are of about 20-30 K, equiva- 
lent to the 2% uncertainty in the IR fluxes of Vega claimed by 
Bohlin & Gilliland (2004) to be a realistic value. 

6. Conclusions 

We have presented a method (called SEDF) to compute ef- 
fective temperatures, angular semi-diameters and bolometric 
corrections from 2MASS photometry. We have adopted an ap- 
proach based on the fit of the observed VJHK magnitudes us- 
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Table 2. Bolometric correction as a function of effective temperature for log g =4.5 and different metallicities (valid for dwarf and subdwarf 
stars). Tables for log g =3.5 and log g =4.0 are available in electronic form. 
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-0.285 


1.813 
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1.003 
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1.022 


-0.199 


0.986 


-0.272 


0.957 
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— 


— 


— 


— 
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— 


— 


— 


— 
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0.931 


-0.031 


0.913 


-0.101 


0.883 


— 


— 


— 


— 


6900 


0.019 


0.884 


-0.026 


0.866 


-0.096 


0.838 










7000 


0.023 


0.837 


-0.023 


0.821 


-0.092 


0.793 










7100 


0.026 


0.792 


-0.019 


0.775 


-0.087 


0.748 










7200 


0.029 


0.746 


-0.016 


0.730 


-0.083 


0.704 










7300 


0.032 


0.701 


-0.012 


0.686 


-0.079 


0.660 










7400 


0.034 


0.656 


-0.009 


0.642 


-0.076 


0.617 










7500 


0.037 


0.612 


-0.007 


0.598 


-0.072 


0.574 










7600 


0.039 


0.568 


-0.004 


0.555 


-0.069 


0.532 










7700 


0.041 


0.525 


-0.002 


0.512 


-0.066 


0.490 










7800 


0.043 


0.481 


0.000 


0.469 














7900 


0.044 


0.438 


0.002 


0.427 














8000 


0.046 


0.396 


0.004 


0.385 















ing synthetic photometry, and it yields accuracies around 1% 
in T e ff, 2% in 8, and 0.05 mag in BC, in the temperature range 
4000-8000 K. A zero point offset was added to the synthetic 
photometry computed from the Kurucz atmosphere models to 
tie in our temperature scale with the Sun's temperature through 
a sample of solar analogues. From the application to a large 
sample of FGK Hipparcos dwarfs and subdwarfs we provide 
parametric calibrations for both effective temperature and bolo- 
metric correction as a function of (V - K)q, [m/H] and logg. 
Note that the method presented here has been selected as one of 
the main sources of effective temperatures to characterize the 
primary and secondary targets of the COROT space mission 
(Baglin et al. 2000). Also, it is being currently implemented 



as one of the tools offered by the Spanish Virtual Observatory 
(Solano et al. 2005). 

The resulting temperatures have been compared with sev- 
eral photometric and spectroscopic determinations. Although 
we obtained remarkably good agreement, slight systematic dif- 
ferences with other semi-empirical methods, such as the IRFM, 
are present. This is probably due to the uncertainties in the ab- 
solute flux calibration used by different techniques. It is possi- 
ble that, in spite of the great effort carried out by Cohen et al. 
(2003a) and others to construct a consistent absolute flux cal- 
ibration in both the optical and the IR regions, some prob- 
lems still remain, which introduce small systematic effects in 
the temperatures. However, these effects seem to be as small 
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Fig. 12. BC K - (V - K) fits for four groups of stars with different 
metallicities. The empirical relationships correspond to log g =4.5 and 
[m/H] =-2.0, -1.0, -0.25 and +0.25. 
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Fig. 13. BC K - (V - K) a relationships for log# =4.5 and four different 
metallicities. 



as 20-30 K and could be explained through uncertainties in 
the IR fluxes of about 2%. In conclusion, the results presented 
here strongly suggest that, given the small differences found 
between methods, the effective temperature scale of FGK stars 
(4000-8000 K) is currently established with a net accuracy bet- 
ter than 0.5-1.0%. 
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Table 3. Effective temperatures, angular semi-diameters, radii and bolometric corrections in both V and K bands, with their uncertainties, for 
10999 stars as resulting from the algorithm presented in this paper. 
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